WHITENING OF THE QUARK-GLUON PLASMA 
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Parton-parton collisions do not neutralize local color charges in the quark-gluon plasma as they 
only redistribute the charges among momentum modes. We discuss color diffusion and color con- 
ductivity as the processes responsible for the neutralization of the plasma. For this purpose, we 
first compute the conductivity and diffusion coefficients in the plasma that is significantly colorful. 
Then, the time evolution of the color density due to the conductivity and diffusion is studied. The 
conductivity is shown to be much more efficient than the diffusion in neutralizing the plasma at the 
^SJ , scale longer than the screening length. Estimates of the characteristic time scales, which are based 

on close to global equilibrium computations, suggest that first the plasma becomes white and then 
' the momentum degrees of freedom thermalize. 
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I. INTRODUCTION 



Production of the quark-gluon plasma is expected at the early stage of high-energy nucleus-nucleus collision when 
the energy density is sufSciently high. The experimental data on the so-called elliptic flow [1], which have been 
r^' obtained at the Relativistic Heavy-Ion Collider (RHIC) in Brookhaven National Laboratory, suggest a surprisingly 
Oh, short, below 1 fm/c [2], equilibration time of the system. Understanding of the thermalization process is thus a key 
^ ' issue of the quark-gluon plasma physics. 

In our previous study [3], we have analyzed the local equilibrium of the plasma, which is defined as a state of 
maximal local entropy. Using the kinetic equations with the collision terms of the Waldmann-Snider form, we have 
proved that such a state is generically colorful, i.e. the color four-current is non- vanishing. Thus, the collisions, which 
are responsible for equilibration of the parton momenta, do not neutralize the local color density. Since the color 
d I current is (covariantly) conserved in every collision process, the inter-parton collisions redistribute the color charges 
among various momentum modes but they do not change a local macroscopic color charge. Consequently, if the color 
charges are not homogeneously distributed in the process of the plasma production due to, say, statistical fluctuations, 
the inter-parton collisions will not neutralize the system. On the other hand, the global equilibrium of the quark-gluon 
plasma is locally colorless because of the maximum entropy principle. We assume here that the system does not carry 
a global color charge and that it does not experience an external chromodynamic field. Once the inter-parton collisions 
are not responsible for the neutralization of the local charges, one has to invoke other collective mechanism to whiten 
the plasma. This is the subject of this article. 

Local charges are neutralized due to the currents that flow in the system. We consider the diffusive currents 
generated by the charge density gradient (Pick's law) and the ohmic currents caused by the chromoelectric field 
(Ohm's law) which is, in turn, induced by the charge density. The color conductivity of the quark-gluon plasma has 
been studied for long time [4-11], but only recently the problem has been well understood [12-19]. As far as we 
know, the color diffusion was only briefly discussed in [11]. In all these papers, the plasma near the colorless global 
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equilibrium was studied. We are, however, interested in the plasma that is locally colorful. Thus, in Sec. II we derive 
the diffusion and conductivity coefficients in such a plasma, and then, in Sec. Ill the temporal evolution of the color 
charge density is considered. The ohmic currents are shown to be much more efficient than the diffusive ones in 
neutralizing the local charges. 

If not stated otherwise, we follow here the same conventions and notations as in [3]. 



II. COLOR DIFFUSION AND CONDUCTIVITY COEFFICIENTS 



In this section we derive, using transport theory, the diffusive and ohmic currents in a plasma that is locally colorful. 
The transport equations of quarks, antiquarks and gluons, which form the basis of our analysis, read 



{D° + v D)Q(p, x) - |{E + V X B, Vpg(p, x)} = C[Q, Q, G] , 
(£>o + V D)Q(p, a;) + |{E + V X B, VpQ(p, x)} = C[Q, Q, G] , 
(D° + V • D)G(p, :r) - |{E + V X B, VpG{p, x)} = Cg[Q, Q, G] . 



(la) 
(lb) 
(Ic) 



Q<^i(p,a;) 


= exp 


Q^^(p,x) 


= exp 




= exp 



The (anti-)quark on-mass-shell distribution functions Q{p,x) and Q{p,x), which are Nc x Nc hermitean matrices, 
belong to the fundamental representation of the SU(A^c) group, while the gluon distribution function G(p, x), which is 
a (iV^ — 1) X (N^ — 1) matrix, belongs to the adjoint representation. The covariant derivative = — ig[A^^{x), • • • ], 
the chromoelectric and chromomagnetic fields, E and B, which enter the transport equations also belong to either 
the fundamental or adjoint representation, correspondingly. To simplify the notation, and differently than in [3], we 
use the same symbols D^, D, E and B to denote a given quantity in the fundamental or adjoint representation. 
X = {t, x) denotes the four-position while p is the three-momentum. Because the partons are assumed to be massless, 
the velocity v equals p/|p|- The collision terms C, C and Cg will be discussed later on. 

We are interested in a state close to the colorful local equilibrium. When the effects of quantum statistics are 
neglected, the (on-mass-shell) local equilibrium distribution functions read [3] 

(2a) 
(2b) 
(2c) 

where = {Ep,p), and Ep = \p\ for massless quarks and antiquarks, and for gluons; /3 = 1/T, w'^, /Zb denote, 
respectively, the inverse temperature, hydrodynamic velocity and baryon chemical potential; the colored chemical 
potentials of quarks (jl) and of gluons (jig) obey the relationship 

Jlg{x) = 2T«IV[r»M(a;)] = M„(a;)T« , (3) 

where r", with a = 1, ... , A^^ — 1 are the SU(A^c) group generators in the fundamental and adjoint representations. 



normalized as Tr[T"T^] = ^S"'' and Tr[r"T''] = A^^'^"'' 



In the local equilibrium state there is a non-vanishing color charge density, which is defined as 



P{x) = -^ J ^{q{p,x)-q{p,x) + 2T'^Tr[T'^Gip,: 



(4) 



with 



q{p, x) = Q{p, x) - -^Tr [Q{p, x)] , q{p, x) = Q{p, x) - -^Tr[Q{p, x)] . 

In the local rest frame, where = (1, 0, 0, 0), the color current defined as 

d^p 



j(a;) = -f / ^v(g(p,x)-g-(p,ar) + 2r«Tr[TO(p,ar) 



(5) 
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vanishes because the distribution functions are locally isotropic. 

We now study the system for long time scales. We consider small deviations from local equilibrium, and we write 
down the quark distribution function as Q{p,x) = Q'^'^{p,x) + 6Q{p,x). Assuming that 

> \5Q\ , \D°Q'^\ > \D°SQ\ , {TtQ^'^l > \T>5Q\ , IVpQ'^^l > \VpSQ\ , (6) 

and taking into account the local isotropy of the equilibrium state, the transport equation (la) can be approximated 
as 

(2?° + V • D) + |{E, VpQ«i} = L[5Q] , (7) 

where L[SQ] is the collision term linearized around the local equilibrium distribution function. Analogous equa- 
tions hold for antiquarks and gluons. Let us recall here that the collision terms evaluated for the local equilibrium 
distribution functions (2) vanish [3]. 

To get the transport coefficients of color diffusion and conductivity we assume that the linearized collision terms 
L, L and Lg satisfy the relationship 

/|f3Vira = -7/|^v*<3, (8) 

where I/7 is the characteristic relaxation time which, for simplicity, is assumed to be the same for quarks, antiquarks 
and gluons. Our analysis based on the linearized transport equation (7) with Ansatz (8) is valid for t > I/7. The 
relation (8) is obeyed [16] by the Waldmann- Snider collision term linearized around global (colorless) equilibrium. 
Such a linear collision term turns out to be non-local in velocities, allowing for the covariant color current conservation 
[12]. A similar linearization around the local colorful equilibrium is, in principle, feasible, but it appears rather difficult 
as it requires knowledge of the scattering matrix elements computed in the colorful background. 
The relationship (8) is trivially satisfied by the collision term in the relaxation time approximation 

C^TA[g, 0, G] = i {Q'^'^ip, x) - Q(p, x))=~ 6Q{p, x) = L^^^[6Q] , (9) 

with r = 1/7. Unfortunately, this approximation is known to contradict the covariant current conservation as C, C 
and Cg in the form (9) violate the collisional invariant 

/ (0^ " ^ + 2r»Tr[T»C,]) = . (10) 

However, the collision term (9) can be improved to a form that complies with the condition (10). In analogy to the 
so-called BGK collision term [20], we have found the expression 

C^GK[Q,Q,G] = ^(Qip,x)-N{x) (7V^<i(a;))-'Q«i(p,ar)) (11) 
« -l.(^SQ{p,x) - 6N{x) {N-'i{x)y'Q''i{p,x)) = L''^'^\bQ\ , 

with 

N{x)^ J ^Q{p,x), (12) 

and (A^°i)^^ being the inverse matrix of N^"^. One easily shows that the collision terms of the form (11) satisfy the 
constraint (10), as 

f d^P ^BGK ^ f d^P ^BGK ^ f d^P ^BGK ^ r. 
J (277)3 J (27r)3 J (27r)3 

The collision term (11) also obeys the relation (8). It still violates the energy-momentum conservation law, but it can 
be further improved [20]. 

Using the relation (8), the color current generated by deviations from equilibrium is 
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iix) = -fl ^v(5g(p,ar) -5g-(p,ar) +2T«TV[T«5G(p,ar))]) (13) 

One observes that the term with drops out as the color current vanishes in local equilibrium. We also take into 
account that the local equilibrium is isotropic and we perform partial integrations in the terms with the chromoelectric 
field. Putting additionally = 1, the current gets the form 

j= 1^^/ (03(9'^'^-5^'^ + 2r"Tr[™^^]) (14) 



.9^ a f d^P 1 
^ 37 y (27r)3 Ej, 

j = _pDp + i (K, E} - ^Tv[a, E]) + T«Tr[T«{a„ E}] , (15) 



which can be written as 
where the diffusion constant V and conductivity coefficients cr^, ag are 



(16) 



As seen, the current (15) is traceless as it should. It gets a much simpler form in the adjoint representation. Namely, 
for j" = 2Tr[T°j], we get 

j"(a;) = -X'D«V^(a;) + a'^'^E^x) , (17) 
where D"*" = ^''''V + gf^'^^'A", and the color conductivity tensor reads 

^''' = ^1 (Tr[{T", r'}{Q''^ + Q'^i)] + Tr[{T«, T^jG''^]) . (18) 

The diffusion constant is, as previously, 1/87. When the equilibrium is colorless, the conductivity a is proportional 
to the unit matrix in color space, but for a colorful configuration it is not. 

It is interesting to note that the chromoelectric field contributes to the induced baryon current. Repeating the 
analysis fully analogous to that of the color current, one finds that the baryon current defined as 

= 3 j ■ 

equals 



b(a;) = -VVb{x) - ^Ti[ag E{x)] , (20) 

where b is the baryon density while T> and (Tg are given by Eqs. (16). When the equilibrium is colorless, the conductivity 
a is proportional to the unit matrix in color space, and the effect of the chromoelectric field on the baryon current 
disappears. 
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III. TEMPORAL EVOLUTION OF THE COLOR DENSITY 



Our aim here is to discuss how a locally colorful quark-gluon plasma becomes white. As an introduction to 
our chromodynamic considerations, we first discuss the temporal evolution of the electric charge density in the 
electromagnetic plasma. 



A. Diffusion vs. conductivity - electrodynamic case 

The electric current (j) generated by both the gradient of charge density (p) and the electric field (E) is 

i{x) = -VVp{x) + aE{x) . 



(21) 



V and a are assumed here to be the transport coefficients derived in a semi-static limit as in Sec. II. Therefore, 

Eq. (21) holds for slowly varying p(x) and E(x). 

Taking into account the Gauss law VE(a:) = p{x), the current conservation dp/dt + Vj = combined with Eq. (21) 
provides the equation 



p{x) = 



which, supplemented by the initial condition p(0,x) = po{x), is solved by 



with n satisfying the diffusion equation 



p{x) 



'■n{x) 



n{x) = , 



(22) 



(23) 



(24) 



and the initial condition n(0,x) = po(x). Eq. (22) can be easily solved by means of the Fourier transformation as 



p{x) 



(27r)3 



-(cr+X>k'')t+ikx 



Po(k) , 



where /9o(k) is the Fourier transform of the initial charge density 

Po(k) = j 



d^xe-*"po(x) . 



(25) 



(26) 



It is assumed here that the integral (26) exsists which requires vanishing of /9o(x) when |x| 
We note that the solution (25) can be also written down as 



oo. 



p{x) = 



1 



d^x' exp { -at- ^^^^^ ) Po(x') 



where the Green's function 



G(x,x',t) 



(47rX>t)3/- 



rexp 



(x - xQ ^ 
4Vt 



at 



(27) 



(28) 



represents the charge density which obeys Eq. (22) and equals d''^\'x. — x') at t = 0. 

The solution (25) shows that the charge density modes of all k decay exponentially. The long-wavelength modes 
with k^ < a/D arc dominantly neutralized by the ohmic currents while those with k^ > a/D are neutralized due to 
the diffusion. It should be remembered, however, we can trust the solution (25) or (27) only for sufficiently long time 
intervals because Eq. (21) holds for slowly varying p{x) and E(a;). We return to this point in the next section where 
it is discussed quantitatively in the context of quark-gluon plasma. 
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B. Diffusion vs. conductivity - chromodynamic case 



The color current generated by both the gradient of color density and the chromoelectric field is given by Eq. (17). 
The covariant current conservation Z)°p + Dj = combined with the Gauss law DE(a;) = p{x), provides the equation 

- X» + 0-] p{x) = , (29) 

where the term [D, cr] E has been neglected. In the Appendix we show that [D, a] E <C ctDE in the small coupling 
limit. All chromodynamic quantities discussed in this section belong to the adjoint representation, and thus, the color 
indices are suppressed. 

Eq. (29) can be treated as its Abelian counterpart (29) if [-D°,cr] = 0. In the Appendix this commutator is shown 
to be indeed small. Thus, Eq. (29) is solved by 

p{x) = e-"* n{x) (30) 

with n satisfying the diffusion equation 

[D° -VB^]n{x) = . (31) 

p{x) and n{x) obey the initial condition p(0,x) = n(0,x) = po{x.). Thus, we expect that, as in the electromagnetic 
case, the charge density decays exponentially and the conductivity dominates over the diffusion for the modes with 
k2 < a/V. 

For further discussion one needs an estimate of the relaxation time I/7 which controls both <j and "D. However, 
the reliable estimate can be given only for the quark-gluon plasma close to global (colorless) equilibrium of very high 
temperatures where 1/g ^ 1. Then, the color conductivity is of order [9,16] 

According to Eq. (18), the conductivity, due to the dimensional argument, can be approximated as cr ~ g^T"^/^ which 
combined with the estimate (32) provides 

^~5Ml/5)T~*^°f*' 

where tgoft is the characteristic time scale of the of parton-parton collisions at momentum transfers of order g'^T [16]. 
We also observe that 

^ ~ g^T^ ~ ml , (34) 

where m^i is the screening mass. Having these estimates, we first note that all modes of charge density longer 
than the screening length are neutralized dominantly by the ohmic currents. However, we can trust Eqs. (30,31) 
only for time intervals longer than I/7 because the derivation of the color conductivity presented in Sec. II is valid 
for t > 1/7. Taking into account the estimate (32), we find that all modes of charge density vanish at i > I/7. 
Since the characteristic time scale of color dissipation cannot be shorter than tgoit, the whitening of the quark-gluon 
plasma occurs at t ~ igoft- At shorter times scales the color density is expected to oscillate. We note that in the 
electromagnetic plasma local charges are neutralized very fast, but the currents survive in the system for a long time 
as the plasma is a very good conductor. Although, the quark-gluon plasma is a rather poor color conductor, the color 
currents can still persist in the system significantly longer than the color charge density [12], as they couple to the 
non-perturbative chromomagnetic fields. 

Since the conductivity is responsible for whitening of the quark-gluon plasma in the long-wave limit, we discuss in 
more detail the equation 

[D°+a]p{x)=0, (35) 
which describes how the ohmic currents neutralize the system. Eq. (35) is solved by 

p{x) = fl{x, xo) e~'^Vo(x) 0(a;o, a;) , (36) 
where x = (f, x), xo = (0,x) and Q,{x,xo) is the parallel transporter 
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t 

n{x,xo) =Texp ig I dt'A'^{t',x) , (37) 
with T denoting the time ordering. Observing that 

- igA^ix)) n{x, xo) = n{xo, x) ( ^ +i<?^°(a;)) = , (38) 

one shows by direct calculation that the expression (36) solves Eq. (35). Since a has non-diagonal entries, various 
colors are coupled to each other in the course of temporal evolution. 

IV. DISCUSSION 

As discussed in our previous paper [3], parton-parton collisions thermalize the momentum degrees of freedom but 
they do not neutralize the local color charges. To whiten the quark-gluon plasma collective phenomena are required. 
T he lo cal charges generate chromoelectric fields, which, in turn, induce color currents. At the scale longer than 
^/V/a, these ohmic currents effectively neutralize the system, more effectively than the diffusive currents caused by 
the charge density gradients. 

The question arises what is the characteristic time scale of momentiim thermalization and that of the plasma 
whitening. Our analysis implicitly assumes that the equilibration of momentum is much faster than the neutral- 
ization, as we linearize the transport equations around the local equilibrium distribution functions which cancel the 
collision terms. In other words, it is implicitly assumed that the plasma momentum distribution first reaches its local 
equilibrium form, and then the system is neutralized. Unfortunately, we are unable to compute the two scales of 
interest as it requires an analysis of parton-parton scattering in a colorful non-equilibrium configuration. Our choice 
of local equilibrium configuration found in [3] is to some extent dictated by technical reasons. The local equilibrium 
distribution functions represent a nontrivial colorful configuration that is convenient to compute transport coefficients 
as the collision terms then vanish. 

The problem of plasma equilibration is also complicated by the fact that color collective phenomena are not only 
responsible for the whitening but they also contribute to the momentum equilibration. One of us has argued for 
long time [21 24], see also [25,26], that color plasma instabilities, which occur in anisotropic systems, speed up the 
momentum thermalization. If the plasma momentum distribution is strongly elongated in one direction, as it occurs 
in heavy-ion collisions, the instabilities generate momentum in the transverse direction, making the system more 
isotropic. However, the instabilities also generate local color charges that have to be neutralized. Thus, the whole 
process of equilibration of the quark-gluon plasma is very complex, and it depends on the plasma initial state. 

The only reliable estimates of the time scales of interest have been found for the perturbative quark-gluon plasma 
which is close to global equilibrium. To equilibrate the system's momentum degrees of freedom, the parton-parton 
interactions with momentum transfers of order T are needed. Such a transfer can be achieved in a single parton-parton 
collision or as a cumulative effect of many soft scatterings. The time scale of such processes is [16] 

~ ,^ln(V,)T • ^''^ 

As argued in the previous section, the whitening of the quark-gluon plasma occurs at f > isoft- Thus, the plasma 
becomes white first and then the momentum degrees of freedom thermalize as tgoft ^ ^hard- Analogous analysis for a 
colorful background should include the effect of the colored chemical potentials that might alter the above picture. 

At the end, let us recapitulate our considerations. Within the QCD transport theory we have found the conductivity 
and diffusion coefficients in the colorful equilibrium configuration. While the diffusion constant is proportional to the 
unit matrix in color space, the conductivity coefficient has a nontrivial tensorial structure. The macroscopic equation 
describing the temporal evolution of the color charge density has been derived. Its solution shows that the ohmic 
currents dominate whitening of the quark-gluon plasma at sufficiently long scale. 
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APPENDIX A 



In this Appendix we argue that the commutators [D", cr] and [D, a] are small in a perturbative regime, but we first 
show that [Z?^, a] ~ (£'^cr) where is the covariant derivative of rank 2. 
We are interested in the following expression 

D„„, (a^'^Efc) = V + gr'^' A^) (a^'-Eb) . (Al) 

Because of the antisymmetry of the structure constants /"''^, Eq. (Al) can be rewritten as 

D„„-(a"'''E6) = D„„-(a"'''E6) + 3/"^"' AV^^Eb + 5/"^"' A V^^Eb- = (Dg.^/a"'"') E;, + a'^^'Tiw^b' , (A2) 

where 

(D)^J = VS'^'S^'^ + gf'""'5'"^ A" + gf'"^5'"' A" . (A3) 
In matrix notation Eq. (A2) gets the form 

D(c7E) = (7 DE + (Dcj) E , (A4) 

and thus 

[D'',a] = {b''a) . (A5) 

Wc are now going to show that [£''*, cr] is negligible when 1/g ^ 1. We actually demonstrate that (D^o-g), computed 
in the fundamental representation, is suppressed by powers of g when compared to CTqD'^. The same analysis can be 
done in the adjoint representation, and for the gluons, but it requires tedious manipulations with color indices. 

We first compute Do-g. Because of the local isotropy of the equilibrium state, we have 

= 1^ / 1? 5 " = f / 5'^ ■ + • 

where we have used the fact that = 1. Using the transport equation (7), Eq. (A6) is rewritten as 

= ^ / (0:^ ( ~ + ^"'^ + f ~ ^'"^^ + ^^^^^^ + ^['^'^^0 • 

The first term in the r.h.s. of Eq. (A6) vanishes because of local isotropy. The remaining two terms are nonzero but 
DfTgE is seen to be smaller than cTgDE by at least two powers of g. The terms with the collision terms are even more 
suppressed. 

Let us now discuss D^Uq. Using the transport equation (7), one gets 

D'cjq = |! y ^ -1 ( - vV(Q'='^ + Q<=<i) + |{E, V,{Q^^ - Q^^)} + {L[5Q] + L[5Q\)) . (A8) 

The first and the second term in the r.h.s. of Eq. (AS) both vanish because of local isotropy of the equilibrium 
momentum distribution. The third term does not vanish but {D°aq)p is suppressed with respect to (JqD^p by powers 
of g hidden in L and L. 



[1] K. H. Ackermann et al. [STAR Collaboration], Phys. Rev. Lett. 86, 402 (2001) [arXiv:nucl-ex/00090111. 

[2] U. W. Heinz and R F. Kolb, Nucl. Phys. A 702, 269 (2002) [arXiv:hep-ph/01 11075]. 

[3] C. Manuel and St. Mrowczynski, Phys. Rev. D 68, 094010 (2003) [arXiv:hep-ph/0306209]. 

[4] U. W. Heinz, Annals Phys. 168, 148 (1986). 

[5] W. Czyz and W. Florkowski, Acta Phys. Polon. B 17, 819 (1986). 



8 



[6] A. Dyrck and W. Florkowski, Phys. Rev. D 36, 2172 (1987). 

[7] St. Mrowczyiiski, Acta Phys. Polon. B 19 (1988) 91. 

[8] St. Mrowczyiiski, Adv. Ser. Direct. Higli Energy Phys. 6, 185 (1990). 

[9] A. Selikhov and M. Gyulassy, Phys. Lett. B 316, 373 (1993) [arXiv:nucl-th/9307007]. 
[10] A. V. Selikhov and M. Gyulassy, Phys. Rev. C 49 (1994) 1726. 
[11] H. Hciselberg, Phys. Rev. Lett. 72, 3013 (1994) [arXiv:hep-ph/9401317]. 
[12] D. Bodeker, Phys. Lett. B 426, 351 (1998) [arXiv:hep-ph/9801430]. 
[13] D. Bodeker, Nucl. Phys. B 559, 502 (1999) [arXiv:hep-ph/9905239]. 
[14] D. F. Litim and C. Manuel, Phys. Rev. Lett. 82, 4981 (1999) [arXiv:hep-ph/9902430]. 
[15] D. F. Litim and C. Manuel, Nucl. Phys. B 562, 237 (1999) [arXiv:hep-ph/9906210]. 
[16] P. Arnold, D. T. Son and L. G. Yaffe, Phys. Rev. D 59, 105020 (1999) [arXiv:hcp-ph/9810216]. 
[17] J. P. Blaizot and E. lancu, Nucl. Phys. B 557, 183 (1999) [arXiv:hep-ph/9903389]. 
[18] P. Arnold and L. G. Yaflfe, Phys. Rev. D 62, 125014 (2000) [arXiv:hep-ph/9912306]. 

[19] J. M. Martinez Resco and M. A. Valle Basagoiti, Phys. Rev. D 63, 056008 (2001) [arXiv:hep-ph/0009331]. 
[20] P.L. Bhatnagar, E.P. Gross, and M.Krook, Phys. Rev. 94, 511 (1954). 
[21] St. Mrowczyiiski, Piiys. Lett. B 314, 118 (1993). 
[22] St. Mrowczyiiski, Phys. Rev. C 49, 2191 (1994). 
[23] St. Mrowczyiiski, Phys. Lett. B 393, 26 (1997). 

[24] J. Randrup and St. Mrowczynski, Phys. Rev. C 68, 034909 (2003) [arXiv:nucl-th/0303021]. 
[25] P. Romatschke and M. Strickland, Phys. Rev. D 68, 036004 (2003) [arXiv:hep-ph/0304092]. 
[26] P. Arnold, J. Lenaghan and G. D. Moore, JHEP 0308, 002 (2003) [arXiv:hep-ph/0307325]. 



9 



